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Abstract 

We derive an exact expression for the transmission amplitude of a 
particle moving through a harmonically driven delta-function potential 
by using the method of continued-fractions within the framework of Flo- 
quet theory. We prove that the transmission through this potential as a 
function of the incident energy presents at most two real zeros, that its 
poles occur at energies nhui + e* (0 < Re(e*) < fuS), and that the poles 
and zeros in the transmission amplitude come in pairs with the distance 
between the zeros and the poles (and their residue) decreasing with in- 
creasing energy of the incident particle. We also show the existence of 
non-resonant "bands" in the transmission amplitude as a function of the 
strength of the potential and the driving frequency. 

1 Introduction 

Time-dependent potentials in mesoscopic systems have been studied for a num- 
ber of years in connection with electron-phonon interactions El , quantum tun- 
neling time ||, ||, ||, ionization [[| electronic transmission H, |L M, [l0| 
and also in the field of quantum chaos . One of the interesting features of 
localized time-periodic potentials is the presence of resonances or quasi-bound 
" states" , which could be thought of as electrons dynamically trapped by the os- 
cillating potential. This is also a feature common to all multi-channel quantum 
scattering problems . 

The solution of any three-term recursion relation, either in functions or in 
operators, is a continued fraction (C.F.). Such is the case for tight-binding 
hamiltonians |l4|, |l5| |, for the time-independent Schrodinger equation (in dis- 
cretised space) ]16(|, and for harmonic time-dependent potentials, such as an 
atom in a standing-wave laser field |l~7], [l8| or for tunneling in the presence of 
phonons |l| . 



We first prove here, starting from Schrodinger's equation, that the transmis- 
sion amplitude has the structure of a C.F. of functions of the incident energy 
and the strength of the delta potential. This is a considerable advantage over 
the numerical computation of the transmission done before using this kind of 
potential (7[ ||. Our expression allows us to study with greater detail, both 
analytically and numerically, several different features of the transmission that 
had not been noticed or explained before, such as the location of the zero-pole 
resonances of the transmission and the almost periodic behavior of their posi- 
tion, the existence of non-resonant "bands", the dependence of the pole residues 
on energy and the existence of the so-called (in the language of nuclear physics) 
"threshold anomalies" in the transmission. We believe our work gives some 
insight and clarifies issues in the general problem of scattering through harmon- 
ically driven localized potentials. One such potential is the Landau-Biittiker 
potential, for which numerical and analytical studies of the transmission have 
been done j^, [To), showing many similarities with the transmission properties 
of the potential studied in this paper. 

In Section 2 we use Floquet's theorem to derive the equations that couple 
different components of the wave function in a plain wave basis, and then we 
use these equations to find the S-matrix for this potential. In Section 3 we 
solve the equations derived in Section 2 and find the exact C.F. expression for 
the transmission amplitudes. In Section 4 and 5 we study analytically and 
numerically the zeros and poles of the transmission. In Section 6 we briefly 
discuss the "threshold anomalies" in the transmission amplitude. 



2 Scattering Matrix and Transmission Ampli- 
tudes 

A formal treatment of the problem of scattering by a time-periodic potential 
can be found in || Also, a discussion about the sub-space of the Hilbert 
space suitable for the treatment of time-periodic potentials can be found in |2(| . 
The Hamiltonian we consider is 

H( X> t) = -^-^ + V8(x)co S (cjt), (1) 

where [i is the mass of the particle. Even though energy is not conserved, the 
Floquet energy, e, is conserved for this system and takes on a continuous range 
of values in the interval 0<e<?ia;. The Floquet eigenstate with Floquet energy, 
e, takes the form 

oo 

* e (x,t)= iWe" |(£+nM ' (2) 

n— — oo 

Since the potential is zero everywhere except at x — 0, we assume ^ n {x) to be 
of the form 

ip n L {x) = -=(a n e ik " x + d n e- iknX ) for x < 
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^ n R (x) = -^=(c n e ik " x +b n e- ik " x ) fon>0. (3) 

The factor ^= has been included to ensure unitarity of the S-matrix and the 
wave- vectors k n are defined by 



k n = \j^{e + nhuj) (4) 

In this paper n will always have the range — oo < n < oo. The square root 
function has its branch cut on the real energy axis (so that a real energy gives 
a real momentum) and for the energy on the negative real axis we will use the 
Riemann sheet that has Im(k n ) > 0, the so called "physical sheet". On this 
sheet the momentum k n is on the positive imaginary axis for any n < 0. This 
is required to allow for evanescent modes (exponentially decaying on both sides 
of the potential) since they can contribute significantly to the wave-function in 
the neighborhood of the delta function. These evanescent states are also related 
to resonances or quasi-bound states which are known to exist in multi-channel 
problems. To avoid any unphysical exponentially growing states when dealing 
with positive pure imaginary momentum states, we will require that a n = and 
b n = for n< — 1. To study the resonances in the transmission we will allow 
the energy to take complex values and will not necessarily stay on the physical 
sheet. This will be discussed in section 5. 

The Floquet eigenstate, & e (x, t), must be continuous at x — 0. This leads 
to the condition that 

a n + d n = c„ + b n . (5) 

Because of the delta function in the Hamiltonian, the slope of ^ e (x,t) is dis- 
continuous and satisfies the condition 



dx x —q+ dx 
This leads to the condition 



^cos(w£)* £ (0,i). (6) 



x=0~ 



c n + d n - b„ - a n = -2i [h n -x(a n -i + d n -i) + h n (a n +i + d n +i)], (7) 

where 

h n = 2 ^ V (8) 
2h yjk n k n+ i 

We can now combine Eqs. (||) and (^), to obtain the following relations between 
coefficients 

c„ + ih n c n+ i + i/i n _ic„_i = o n - ihnbn+i - ih n -ib n -i (9) 

and 

d n + ih n d n +i + ih n -xd n -\ = b n - ih n a n+1 - ih n -ia n -i- (10) 
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It is useful to separate the propagating modes from the evanescent modes. 
Let us define column vectors 



/c \ 

Cl 

V ; J 



and 



V : J 



(ii) 



where c p contains the amplitudes of the propagating modes outgoing to the 
right, and c e contains the amplitudes of the evanescent modes that decay to 
the right. Analogous definitions apply for the column vectors a p , b p , d p , and 
d e . Note that a e =0 and b e =0. We can now rewrite Eqs. (|9|) and ( [l0| ) in the 
following form 



(1 pp ~t~ XppyCp -\- Xp e 'C e — dp Xpp'bp, 

Xep'Cp ~f~ (lee Ag e )'C e — X e p'bp : 

(lpp Xppj'dp -\- Xp e -d e — bp Xpp'Qp, 

^ep'dp -\- (lee ~t~ A ee )-<i e X e p'(lp, 



(12) 



where l pp and 1 £ 



are infinite dimensional unit matrices and matrices X pp , X ee , 



X ep , and X pe have matrix elements 



(App) m _ m ' 
(Ap e ) r 



* = ih-i S m o6 u _i, 



-i), {X ee ) v . v > - 
(X e p*) urn 



= ih-x <5m,o^,-i- (13) 



Note that we have introduced the convention that the indices m = 0, 1, 2, .., oo, 
and v = —1, —2, —3.., — oo, to help separate propagating modes from evanescent 
modes. 

We can now write the scattering matrix for this system. The scattering ma- 
trix, S, connects the incoming propagating modes to the outgoing propagating 
modes 



S- 



f' 



(14) 



where i and t' are matrices of transmission probability amplitudes and f and f 1 
are matrices of reflection probability amplitudes. More specifically, the matrix 
element t m ^ m i = (t) m ,m' is the probability amplitude for the m' mode entering 



from the left to transmit to the right, and 



is probability ampli- 



tude that the m! mode entering from the left and will reflect to the left. Matrix 
elements of r' and F contain reflection and transmission coefficients for modes 
entering from the right. After some algebra one can show that 



t = t' 



where 



(lpp + Y pp ) 



and r = r' = -(lpp + Y pp ) 1 -Y % 



! PPI 



Ape ' (lg 



Ye,. 



L -A R 



(15) 



(16) 



L pe V ee ~ ee I ^*-ep* 

The effects of the evanescent modes are now explicitly included in the scattering 
matrix. 
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For propagating modes entering in the m channel on the left, the total 
probability for transmission to the right is 

oo 

Tm ^ y |^m'.m| ■ (-^^) 

m'=0 

In the next section we focus on the transmission probability, To. By expanding it 
in continued fractions, we can determine all the transmission zeros and complex 
poles of this system. 



3 Continued fractions solution 

Let us now consider the special case of a single propagating mode, entering in 
channel m = from the left and no propagating modes entering from the right. 
In this case a n — 5 n fi and b n = 0. The probability amplitude for the particle 
to be transmitted into the n th channel (propagating or not) on the right is c„. 
For the propagating modes (n > 0), c n = t n ^. Let us now define the following 
quantity 

f n = — — -co < n < oo (18) 

Cn+l 

In terms of /„, Eq.(||) gives, 

1 + ^ + ihn-ifn-i = for n^O (19) 

In 

1 + ^ + ih-xf-x = — when n=0 (20) 
Jo c 



For n > we can write the solution of (|19|) in the form 

1 ih 

lUn Jn+X 



fn = ^-(l + ^) (21) 



or, 



/n = =^-(l+ % ) (22) 

la n , . ,L n+2 



h 2 

j , "n+3 



For n < — 1 we write the solution in the form 

1 + ihn-xfn-l 



fn = , , .u ik \ (23) 



or, 

/.. - m 

1 H — 

h 2 

1 + . K - 2 



h 2 
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From the above expressions for the /„'s (when n = and n = —1) we can 
obtain cq from Eg. (pp|) , 



c o — to.o 



1 



l + ^ +lh _ lf _ l 1 , *§ , ^-t 



/o ' ^ ft? -r h 2_ 



2 



hi ^ , hi 3 



\ h\ ^ h\ 

'(25) 

With the solution for the /„ s and cq known, any coefficient c n can be found 
(using Eq.[l8|) in the following way, 

for ri>l 



fn-lfn-2—fo 

and 

Cn = fnfn+l — -f-lC for n< - 1 (26) 

From here, the transmission probability, To, can be written as 



oo ^ 

To = ]TM 2 = IM 2 (i + £ — i 

n=0 ra=1 11„'=1 l/n'-ll 



-) = |<o,o| 2 ^(e,^). (27) 



As it can be seen in Fig. none of the f' n s (for n > 0) appear to have zeros in 
the half-plane Re[e] > 0. This implies that the function S(e, V) has neither zeros 
nor poles in that region; consequently, the zeros and poles in the transmission 
probability are the zeros and poles of |io,o| 2 - From now on we concentrate on 
t ,o only. 



4 Transmission Zeros 

The coefficient cq can be taken to be a continuous function of the incoming 
energy E = e + nhuj, instead of a function of the Floquet energy e. This is so 
because, when the incident energy is in channel m (E = mhu> + s), we need 
to solve Eq.^) with the condition a n = 5 n ,m- The solution is given in terms 
of the coefficient c m which now plays the former role of cq. The C.F. for c m 
is given by Eq.(p5|) with m added to all subscripts. As it can be seen easily, 
Cm(s) = Co(mhu> + e) = cq(E), which means that the general solution of the 
problem, for any incoming energy can be obtained from the C.F. expression for 
Co derived in Eq.(|25|) evaluating it at any energy E. Notice also that cq(E) = 
co(mhuj + e) = t. m ^ rn {e). Because of this, and to be rigorous with the notation, 
we define to(E) = t {mhuj + e) = t m . m (e) — Cq(E). From this we can say that 
tg (E) contains the same information as the whole diagonal of the transmission 
matrix t in Eq.Jl^). 



G 



To study the properties of cq(E) it is convenient to define the following 
quantities (see Eq. ^) 

2 " 

g n (e) = h n (e) 



\/e + n\J e + n + 1 
with dimensionless parameters 

mV 2 e 
a = — 5—, e = — - 

and 

e = — = e + n for nftu < E < (n + l)hu) (28) 
Also, we define the function 

F (e) = F (n + e) = F n (e) = 1 + *^±4 = 1 + (29) 

fo{n + e) 1 gn+ngj 

x .gn+2(e) 



Notice that g n (e) also depends on a. 

Using the definitions given above we can write the coefficient co(e) for the 
range n < e < (n + 1), in the following way 

co(e) = c (n + e) = ^- ^ (30) 

F n (e) + ^ 



9o (e) 



In the above expression all quantities are real (for e real) except for g~i(e) which 
is pure imaginary. This has an important consequence for the number of real 
zeros in Co(e) as we shall see next. 

Let us rewrite the equation above in a slightly different way: 

c Q (n + e) = -— (31) 

^™(, e J H / \ 

1 gn-2(f) 
1+ 

l + »G(e) 
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where the continued fraction G(e) is defined as 



G(e) 



a 



(32) 



9-2 (e 



1 + 



Notice G(e) is a real function for e real and < e < 1. 

If Eq. ([u]) is put in the form of a standard fraction, we get 



1 



for n=Q 



Fo+iG 



and 



c (n + e) 



iP n (g n _ 2 ,g n -3, -~,go)G + R n (g n— 2 ) 9n — 3- • • j 90} 
iQn(F n , <?„_i, ...,go)G + S n {F n , g n -i, ■••,9o) 



for n>0 (33) 



In the last expression, P n , Q n , R n , S n are polynomials on the variables g n and F„ 
indicated in parentheses, such that the coefficient of every term is equal to one 
and all variables appear elevated to the first power only (i.e. 1 + g\ + g\g2 + ••)■ 
Notice that G is the only function that can take values between —00 to 00. The 
other functions (F n and g n ) are strictly positive and finite for e > . This 
implies that P n ,Qn, Rn, S n are strictly positive and finite, which means that, 
for n > 0, the numerator of Co can never be zero and the denominator does not 
go to infinity unless G — > 00, in which case the numerator would also blow up 
keeping the fraction strictly positive. From this we can conclude that there are 
no real transmission zeros for incident energies E > hio. 

For n=0 we see that a real zero can only happen when _fo( e : a ) — y 00 or 
G(e, a) — ► 00. The first case only occurs when e = (G blows up at this point 
too). For the second case, the real zeros of t (e,a) are given by the zeros of 
G _1 (e, a). As it can be seen from Fig. ^ the function G _1 (e, a) seems to have 
some periodicity. This can also be seen in Fig. || where the curves 5(a) that 
satisfy to(S, a) = G^ 1 (5, a) — are shown. This dependence of the transmission 
real zeros with the parameter a can be seen directly in Fig. ^ in the sequence of 
transmission graphs for different values of the strength of the delta. 
It is interesting to notice in Fig. || and Fig. |]that there are intervals of a, around 
integer values, for which the real zero in the transmission disappears. A table 
(for < a < 9) with the exact values (to seven digits) of the intervals for which 
there is a real zero in the transmission (apart from the trivial one at E = 0) is 
shown in Table [l[ This behavior of the real zero seems to be matched by the 
behavior of the poles of the transmission as we will see next. 

5 Transmission Poles 

In this section we prove by induction that the poles of co(e) occur at energies 
e = (n + e*) with < Re(e*) < 1 and show that their residue decreases with 
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increasing n (energy). We start with equation (pCl) and write it in the form 



C0in + e) =F n (e) + G n (e) (34) 



where F n and G n satisfy 



F n (e) = 1 



9n(e) 
F n +l(e) 



^'^TTc^h (35) 

Let's assume that co(e) has a pole at e = n + e*. This implies that cp(e) has 



also a pole at e = n + 1 + e* because we can write from Eqs. ( J34[ ) and ( p5| ) 

c (n + 1 + e) = — TaTt T7\ 

1 



9n{e) flw(e) 



F n (e)-1 1 + G„(e) 

= (l + G n (e))(F n (e)-l) 
S„(e)0F„(e) + G„(e)) 

= (1 + Gn(0) (36) 

F n+1 (e)(F n (e)+G n (e)) 1 J 

It is important to notice that, since /o(e) (as shown in Fig.[l]) has neither poles 
nor zeros in the half plane Re(e) > 0, the function F n (e) = F (n + e) = 1 + i h f °\ e \ 
does not have poles and ^(e) ^ 1. This implies that (from the last line of 
Eq. (|(|) if c (e) has a pole at n + e* (which means that F n (e*) + G n (e*) = 0), 
then it must have a pole at n + 1 + e* (The numerator in Eq. ( |36"| ) does not 
vanish at e = e* because G n (e*) — — F n (e*) ^ —1). 

Notice also that when n — > oo, F n — ► 1, because fo(n + e) — > oo, and 
/io(n + e) in this limit. This means that the location of the transmission 
zero approaches the location of the pole as the incoming energy grows; the zero 
happens when e = S, with G n (S) — — 1; the pole occurs when e = e*, with 
G„(e*) = — F n (e*). Obviously S — > e* because F n (e*) — ► 1 as n — > oo. From 
that we conclude that the residue of the poles in the transmission amplitude 
tend to zero as n — ► oo. This explains why, even though the transmission 
amplitude has an infinite number of poles separated by a distance of hu in the 
incoming energy, all at the same distance from the real axis, only the effect 
of the first poles can be seen in the graphs of transmission probability versus 
incoming energy. Another way to say this is that the zeros (in the complex 
plane) at higher energies are very close to the poles, therefore canceling out the 
possible effect of the poles in the transmission. 
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In Fig. H we show a graph of the imaginary part of t (e) where the poles in 
channels 0, 1, 2, can be seen (we refer to channel n as the strip in the complex 
energy plane that satisfies n < Re{e) < n + 1). It is evident in that graph that 
the poles have support on different sheets; we will call these sheets S n , where n 
refers to the channel number. This sheeted structure comes from the fact that, 
because the functions k n (e) have a branch point at e = — n and two Riemann 
sheets, any function of k n (e) (io(e) in particular) will have a multiple-sheeted 
structure. What is called the "physical" sheet (P) in the context of multiple 
channel scattering is the sheet obtained when selecting all the Riemann sheets 
with Im(fc„) > 0. Each sheet of the full multi-sheeted surface can be labeled by 
the sequence of signs of lm(k (X> ),..., Im(fci), Im(fco), Im(fc_i),..., Irr^fc-oo)., e.g., 
(+.. — (bold has been used to indicate the sign of Im(fco)); fortunately 
we do not need to consider all these sheets, since only a small fraction of them 
are of physical interest; the principal one being P or (+..H — Crossing 
the real axis from P at an energy m < e < m + 1 crosses all the branch cuts 
whose branch point is at an energy smaller than e — m (see figure |6|) . The 
sheet that is the smooth continuation of P into the lower complex plane (smooth 
even at the cut) is the one obtained by taking Im(fc_ n )< for all n < m and 
Im(fc_„)> for all n > m . We call these sheets S m and they are precisely the 
ones where the poles in the transmission are found, as shown schematically in 
Fig. 0. 

Resonances, as opposed to bound state poles, do not occur in the "physical" 
sheet; however, from the unitarity of the S-matrix extended to the complex 
plane, they can have important effects on the behavior of the S-matrix along 
the real axis (see jllj for an excellent discussion on multiple channel scattering). 

In single channel scattering, the unitarity of the S-matrix on the real axis 
has consequences in the analytic structure of the S-matrix in both sheets, the 
most important one being the zero-pole structure of the resonances: every pole 
has a companion zero at a position complex conjugate to the position of the 
pole but on a different sheet. For multiple channel scattering the position of the 
pole and zero are not so simply related and often the pole and zero appear in 
the same sheet. 

The position of the pole (e*) in the first channel, can be found by looking 
at the zeros of the function Fo(e, a) + Go(e,et) with the appropriate selection 
of Riemann sheets, as described before, so as to be on the sheet So. A graph 
of the position of the pole as the strength of the potential is varied is shown in 
Fig. ||, where it can be seen that the pole approaches the real axis at e* = 1 
when a — > 0, annihilating the real zero which also goes to the same position 
with a — > 0. As the parameter a is increased the zero moves to the left on the 
real axis (see Fig. H) and the pole moves away from the real axis describing 
an arch. At a = Oq = 0.7821147 (all values of a are accurate to the seventh 
digit) the zero disappears at 8 = and so does the pole. For values of a slightly 
greater, the transmission has a very high peak but it is not a pole, as can be 
seen in Fig. || (compare with the pole at a = 0.78 in Fig. |l0|). In the interval 
Oq < a < a\ the transmission does not have any zero-pole resonances. In this 
interval the peak or false "pole" follows a trajectory that seems the continuation 
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of the trajectory the pole had followed until it vanished. When a = a\ the zero 
reappears on this channel, entering at 5 = 1. At this value however, we do not 
see any pole in this channel. There is in fact a pole on 5*0 but it is located in 
the next channel (1 < Re(e*) < 2). For a « 1.3 the pole has finally made it 
into the n = channel. As a is further increased the pole continues to describe 
an arch until it disappears along with the zero. A similar behavior to the one 
described above occurs for increasing values of the parameter a, with the pole 
describing ever increasing arches, each one farther away from the real axis than 
the previous one. An identical behavior of the poles on other sheets is expected 
from the discussion following Eq.(|36|). 

This quasi-periodic behavior of the transmission resonances as a function of 
the strength of the oscillating potential and the driving frequency has not been 
found before. Around integer values of a the oscillating potential seems to be 
incapable of dynamically trapping particles not even for a short period of time. 

6 Threshold anomalies 

One of the most evident features in the transmission versus energy curves for 
all values of a is what seems to be a discontinuity in the slope at e = nhu 
(channel openings). A close look into this regions reveals that there is a rapid 
divergence of the slope at the thresholds, that is, ^ — > ±oo as e — ► nhui. 
This is actually a common occurrence in all multi-channel scattering problems, 
and it can be traced back to the fact that when the energy is in the n th channel 
near threshold, the first open-channel momentum is = \J e — nfwj. It can be 
proven that the S-matrix has elements that are linear functions of this momen- 
tum near threshold. This clearly implies that the derivative of these elements 
with respect to the energy must diverge at threshold (for more details see pl[|). 
These threshold anomalies can easily be proven to exist in our particular time 
dependent potential by looking at our C.F. solution for to(e) in Eq. [2^. From it, 
it is clear that the derivative of this expression gives an infinite number of terms, 
each one proportional to the derivative of some function g n (e). These functions 
and their derivatives diverge at their thresholds (see the expression for g n in 
Eq. ^8j), therefore the derivative of t (e) must also diverge at each threshold. 
According to this, the shape of the threshold anomaly can be of four different 
kinds, two cusps like and two rounded steps, as shown in Fig. [n]. These four 
different kinds of anomalies are shown in Fig. [l2] as they occur in the |fo(e)| 
graphs for two different values of a. 

7 Conclusions 

The exact solution for the problem of transmission of a particle through a 
monochromatic oscillating potential has been found with the use of C.F's. This 
approach represents a significant improvement over the standard method which 
can only be carried out numerically since it involves the inversion of an infinite 
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tri-diagonal matrix, which can be a computationally demanding problem when 
studying the zero-pole structure of the transmission on the complex plane. The 
efficiency of the C.F. method allowed us to perform, in a short time, calcula- 
tions that required including processes of phonon emission and absorption of 
65th order(for a = 9) which would have required in the standard method the 
computation of a 131x131 S-matrix for every point in the complex plane where 
the transmission is to be evaluated. Most importantly, it allowed us to prove 
rigorously some general properties of the transmission for this system, such as 
the existence of transmission zeros only in the first channel, the location of 
the poles at regular intervals of fiuj in the incident energy, the decrease in the 
residue of the poles with increasing energy, and also allowed us to understand 
the existence of threshold anomalies at the channel openings. 

The study of the behavior of the poles and real zeros as a function of the pa- 

v 2 

rameter a — showed an interesting quasi-periodic dependence on a and the 
presence of " bands" of non-resonant values of a for which the transmission does 
not have zero-pole resonances. It can also be mentioned that the life-time of the 
quasi-bound states decreases with increasing a, which can be accomplished by 
either increasing the strength V of the potential or by decreasing the frequency 
of the oscillations. 

The scattering of particles through a time-periodic potential is another ex- 
ample of multi-channel scattering, for which a great deal of theory and research 
has been done in the past and still continues to be an important topic of re- 
search. In this paper we have put Floquet scattering in the greater context of 
multi-channel scattering to which it belongs and by doing so we believe some 
aspects of the transmission through oscillating (localized) potentials have been 
clarified. 

Even though these results where derived for an oscillating delta potential, a 
similar behavior can be expected to be found for other oscillating finite-range 
potentials such as the Landauer-Biittiker potential, since this potential reduces 
to the oscillating delta in an appropriate limit, and also because several common 
features can already be observed when comparing the transmission curves for 
this two potentials (see ||) . The exact solution for the oscillating delta plus a 
static (attractive or repulsive) delta potential can be found using the same C.F. 
method. The transmission through this potential has been studied numerically 
before Q], showing features similar to the ones discussed here with the main 
difference been that the location of the resonances is displaced by an amount of 
energy roughly equal to the energy of the bound state associated with the static 
part of the potential. 
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I 


h 


a 





0.7821147 


ai 


1.1652568 


1.7710590 


a 2 


2.1873508 


2.7667368 


A3 


3.1979937 


3.7642963 


ai 


4.2045658 


4.7626808 


a 5 


5.2091415 


5.7615115 


a 6 


6.2125623 


6.7606150 


a 7 


7.2152455 


7.7598995 


a s 


8.2174231 


8.7593115 



Tabic 1: Intervals of &= " l X for which a zero-pole resonance can be found in 
the system. The notation a l n /a^ (used in the text and in some graphs) refers to 
the lower/higher value of the interval a n . 
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lm < e ) Re(e) 



Figure 1: Graph of |/o(-E)| 2 in the complex plane (E in units of hui). a = 10. 
Notice that f n {z)=fo{n + e)=fo(E), therefore this graph shows the behavior of 
/ n (e) for < n < 9. The function behaves similarly, but with a different scale, 
for other values of a. 




Figure 2: Graph of the function \G 1 (e, a)\. The zeros of this function give the 
real zeros of the transmission. 
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Figure 3: Location of the real zero of to as a function of a. 
t (S,a) =G~ 1 (S,a) =0. 
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Figure 4: Sequence of |io(e)| graphs for increasing values of a showing the 
evolution of the transmission real zeros. 
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Figure 5: Graph of the imaginary part of to(e) for a=0.5, showing several poles. 
Each pole is located on a different sheet. 
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Figure 6: top view of the physical sheet P and the threshold branch points 
and cuts in the complex energy plane. Indicated with a thin line are the branch 
cuts corresponding to the branch points n — 0,1,2; the thick line represents 
all the branch cuts with branch point at negative energy(n < branch points). 
Paths that start on the upper half of P and go under the branch cut lead to the 
different sheets S n (which in this figure are assumed to be underneath P). 
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Figure 7: Sketch of Im(t{j(e)) versus Re(e) (we assume Im{e) — constant < 0) 
which gives a schematic representation of the different sheets where the poles 
of the transmission are located. The peak structures in the lower part of the 
figure represent poles. 
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Figure 8: Trajectory of one of the transmission poles as a is changed from to 
9. Here a l n /a^ refers to the lower/higher value of the interval a n given in Table 
1. 
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Figure 9: false "pole" of to(e) for a = 0.79, right after the zero has disappeared 
at S = for a = aft. 
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Figure 10: True pole of the transmission for a = 0.78, slightly lower than the 
value of a (a = aft) for which the zero disappears. 
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Figure 1 1 : Various forms of transmission amplitude behavior at a channel open- 
ing. 
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Figure 12: Graph of |to(e)| vs. incoming energy. Threshold anomalies of "cusp" 
type at the channel openings for a=0.5 . Threshold anomalies of "rounded step" 
type for a=2.05 . 
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